The ChEMBL Database in 2023: a drug discovery platform spanning multiple bioactivity data types and time periods.

ChEMBL (https://www.ebi.ac.uk/chembl/) is a manually curated, high-quality, large-scale, open, FAIR and Global Core Biodata Resource of bioactive molecules with drug-like properties, previously described in the 2012, 2014, 2017 and 2019 Nucleic Acids Research Database Issues. Since its introduction in 2009, ChEMBL's content has changed dramatically in size and diversity of data types. Through incorporation of multiple new datasets from depositors since the 2019 update, ChEMBL now contains slightly more bioactivity data from deposited data vs data extracted from literature. In collaboration with the EUbOPEN consortium, chemical probe data is now regularly deposited into ChEMBL. Release 27 made curated data available for compounds screened for potential anti-SARS-CoV-2 activity from several large-scale drug repurposing screens. In addition, new patent bioactivity data have been added to the latest ChEMBL releases, and various new features have been incorporated, including a Natural Product likeness score, updated flags for Natural Products, a new flag for Chemical Probes, and the initial annotation of the action type for ∼270 000 bioactivity measurements.

Exploring the conformational space of bridge-substituted dithienylcyclopentenes.

Stimuli responsive compounds and materials are of high interest in synthetic chemistry and materials science, with light being the most intriguing stimulus due to the possibility to remote control the physicochemical properties of a molecule or a material. There is a constant quest to design photoswitches with improved switching efficiency and especially diarylethene-type switches promise photo cyclization quantum yields up to unity. However, only limited attention has been paid towards the influence of the solution conformation on the switching efficiency. Here, we describe a detailed NMR spectroscopic investigation on the conformational distribution of bridge-substituted dithienylcyclopentenes in solution. We could discriminate between several photoactive and photoinactive as well as two diastereomorphous conformations and show that the trends observed in the switching efficiency match the conformer populations obtained from state of the art NMR parameters in solution.

A combination of spin diffusion methods for the determination of protein-ligand complex structural ensembles.

Structure-based drug design (SBDD) is a powerful and widely used approach to optimize affinity of drug candidates. With the recently introduced INPHARMA method, the binding mode of small molecules to their protein target can be characterized even if no spectroscopic information about the protein is known. Here, we show that the combination of the spin-diffusion-based NMR methods INPHARMA, trNOE, and STD results in an accurate scoring function for docking modes and therefore determination of protein-ligand complex structures. Applications are shown on the model system protein kinase A and the drug targets glycogen phosphorylase and soluble epoxide hydrolase (sEH). Multiplexing of several ligands improves the reliability of the scoring function further. The new score allows in the case of sEH detecting two binding modes of the ligand in its binding site, which was corroborated by X-ray analysis.

End-organ protection in hypertension by the novel and selective Rho-kinase inhibitor, SAR407899.

AIM: To compare the therapeutic efficacy of SAR407899 with the current standard treatment for hypertension [an angiotensin converting enzyme (ACE)-inhibitor and a calcium channel blocker] and compare the frequency and severity of the hypertension-related end-organ damage. METHODS: Long-term pharmacological characte-rization of SAR407899 has been performed in two animal models of hypertension, of which one is sensitive to ACE-inhibition (LNAME) and the other is insensitive [deoxycorticosterone acetate (DOCA)]. SAR407899 efficiently lowered high blood pressure and significantly reduced late-stage end organ damage as indicated by improved heart, kidney and endothelial function and reduced heart and kidney fibrosis in both models of chronic hypertension. RESULTS: Long term treatment with SAR407899 has been well tolerated and dose-dependently reduced elevated blood pressure in both models with no signs of tachyphylaxia. Blood pressure lowering effects and protective effects on hypertension related end organ damage of SAR407899 were superior to ramipril and amlodipine in the DOCA rat. Typical end-organ damage was significantly reduced in the SAR407899-treated animals. Chronic administration of SAR407899 significantly reduced albuminuria in both models. The beneficial effect of SAR407899 was associated with a reduction in leukocyte/macrophage tissue infiltration. The overall protective effect of SAR407899 was superior or comparable to that of ACE-inhibition or calcium channel blockade. Chronic application of SAR407899 protects against hypertension and hypertension-induced end organ damage, regardless of the pathophysiological mechanism of hypertension. CONCLUSION: Rho-kinases-inhibition by the SAR407899 represents a new therapeutic option for the treatment of hypertension and its complications.

CROSS: an efficient workflow for reaction-driven rescaffolding and side-chain optimization using robust chemical reactions and available reagents.

A novel procedure (CROSS: Computational Rescaffolding and Optimization using Synthetic Schemes) for in silico rescaffolding and side-chain optimization is reported with explicit consideration of the route of synthesis and availability of compatible chemical reagents. We have defined a set of retrosynthetic disconnections representing robust reactions, amenable for combinatorial chemistry. This rule set is used to generate virtual fragment databases from available reagents. Each reactive center is annotated with its compatibility with regard to the chemical reactions. The rule set is then applied to a new molecule to obtain separate query subunits for rescaffolding by 3D similarity searching in specific reagent-derived fragment databases. Thus, only fragments compatible with the chemistry and shape of the corresponding query moiety are investigated further. The identified fragment hits directly indicate (1) available chemical reagents that can replace the query moiety in the starting molecule and (2) the route for the synthesis of the proposed molecules.

Accounting for conformational variability in protein-ligand docking with NMR-guided rescoring.

A key component to success in structure-based drug design is reliable information on protein-ligand interactions. Recent development in NMR techniques has accelerated this process by overcoming some of the limitations of X-ray crystallography and computational protein-ligand docking. In this work we present a new scoring protocol based on NMR-derived interligand INPHARMA NOEs to guide the selection of computationally generated docking modes. We demonstrate the performance in a range of scenarios, encompassing traditionally difficult cases such as docking to homology models and ligand dependent domain rearrangements. Ambiguities associated with sparse experimental information are lifted by searching a consensus solution based on simultaneously fitting multiple ligand pairs. This study provides a previously unexplored integration between molecular modeling and experimental data, in which interligand NOEs represent the key element in the rescoring algorithm. The presented protocol should be widely applicable for protein-ligand docking also in a different context from drug design and highlights the important role of NMR-based approaches to describe intermolecular ligand-receptor interactions.

An NMR-based scoring function improves the accuracy of binding pose predictions by docking by two orders of magnitude.

Low-affinity ligands can be efficiently optimized into high-affinity drug leads by structure based drug design when atomic-resolution structural information on the protein/ligand complexes is available. In this work we show that the use of a few, easily obtainable, experimental restraints improves the accuracy of the docking experiments by two orders of magnitude. The experimental data are measured in nuclear magnetic resonance spectra and consist of protein-mediated NOEs between two competitively binding ligands. The methodology can be widely applied as the data are readily obtained for low-affinity ligands in the presence of non-labelled receptor at low concentration. The experimental inter-ligand NOEs are efficiently used to filter and rank complex model structures that have been pre-selected by docking protocols. This approach dramatically reduces the degeneracy and inaccuracy of the chosen model in docking experiments, is robust with respect to inaccuracy of the structural model used to represent the free receptor and is suitable for high-throughput docking campaigns.

pharmACOphore: multiple flexible ligand alignment based on ant colony optimization.

The flexible superimposition of biologically active ligands is a crucial step in ligand-based drug design. Here we present pharmACOphore, a new approach for pairwise as well as multiple flexible alignment of ligands based on ant colony optimization (ACO; Dorigo, M.; Stützle, T. Ant Colony Optimization; MIT Press: Cambridge, MA, USA, 2004). An empirical scoring function is used, which describes ligand similarity by minimizing the distance of pharmacophoric features. The scoring function was parametrized on pairwise alignments of ligand sets for four proteins from diverse protein families (cyclooxygenase-2, cyclin-dependent kinase 2, factor Xa and peroxisome proliferator-activated receptor γ). The derived parameters were assessed with respect to pose prediction performance on the independent FlexS data set ( Lemmen, C.; Lengauer, T.; Klebe, G. J. Med. Chem. 1998, 41, 4502 - 4520) in exhausting pairwise alignments. Additionally, multiple flexible alignment experiments were carried out for the pharmacologically relevant targets trypsin and poly (ADP-ribose) polymerase (PARP). The results obtained show that the new procedure provides a robust and efficient way for the pairwise as well as multiple flexible alignment of small molecules.

Determination of the interfacial water content in protein-protein complexes from free energy simulations.

The question as to how many tightly or weakly bound water molecules are located in interfaces between protein-protein complex constituents is addressed from a phase equilibrium point of view by developing a theory in the canonical ensemble. A fast method based on free energy simulations is described for computing the number of water molecules in the interface regions. Results are given for 211 interfacial cavities of 26 antigen-antibody complexes for which experimentally determined structures are found in the Protein Data Bank. The accuracy of the method is assessed and the computational water content is compared with experimental data, revealing the amount of water molecules not resolved by experimental approaches.